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Abstract 

Classical black hole analogues (alternatively, the analogue systems) are fluid dynamical analogue 
of general relativistic black holes. Such analogue effects may be observed when acoustic perturba- 
tions (sound waves) propagate through a classical dissipation-less tran-sonic fluid. The acoustic 
horizon, which resembles the actual black hole event horizon in many ways, may be generated 
at the transonic point in the fluid flow. Acoustic horizon emits quasi thermal phonon spectra, 
which is analogous to the actual Hawking radiation, and possesses the temperature referred as 
<T^> ■ the analogue Hawking temperature, or simply, the analogue temperature. 

Transonic accretion onto astrophysical black holes is a very interesting example of classical 
analogue system found naturally in the Universe. An accreting black holes system as a classical 
analogue is unique in the sense that only for such a system, both kind of horizons, the electromag- 
netic and the acoustic (generated due to transonicity of accreting fluid) are simultaneously present 
| in the same system. Hence an accreting astrophysical black hole is the ideal-most candidate to the- 

oretically study and to compare the properties of these two different kind of horizons. Also such 
system is unique in the aspect that general relativistic spherical accretion onto the Schwarzschild 
black hole represents the only classical analogue system found in the nature so far, where the 
analogue Hawking temperature may be higher than the actual Hawking temperature. 

o 
o 



o 



1 Black Holes 

Black holes are vacuum solutions of Einstein's field equations in general relativity. Classically, these 
objects are conceived as singularities in space time, censored from the rest of the Universe by mathe- 



qh matically defined one way surfaces, the event horizons. The space time metric defining the vacuum 
exterior of a classical black hole, and the black hole itself, is characterized by only three parameters, 
the mass of the black hole Mbh, the rotation (spin) J and charge q. For J = q = 0, one obtains 
a Schwarzschild black hole, and for q = one obtains a Kerr black hole. These two kind of black 
holes are important in astrophysics. In astrophysics, black holes are the end point of gravitational 
collapse of massive celestial objects. Astrophysical black holes may be broadly classified into two 
categories, the stellar mass (M bh~ a few M & , where M Q is the mass of the Sun), and super massive 
(M BH >10 6 M e ) black holes (SMBH). While the birth history of the stellar mass black holes is theo- 
retically known with almost absolute certainty (they are the endpoint of the gravitational collapse of 
massive stars), the formation scenario of the supermassive black hole is not unanimously understood. 
A SMBH may form through the monolithic collapse of early proto-spheroid gaseous mass originated 
at the time of galaxy formation. Or a number of stellar/intermediate mass black holes may merge 
to form it. Also the runaway growth of a seed black hole by accretion in a specially favoured high- 
density environment may lead to the formation of SMBH. However, it is yet to be well understood 
exactly which of the above processes routes towards the SMBH formation, see, e.g. Rees 2002 for a 
comprehensive review on the formation and evolution of SMBH. Both kind of astrophysical black 
holes, the stellar mass and SMBH, however, accrete matter from the surrounding. Depending on 
the intrinsic angular momentum content of accreting material, either spherically symmetric (for zero 
angular momentum flow), or axisymmetric (for flow with non-zero finite angular momentum) flow 
geometry is invoked to study an accreting black hole system (Frank, King & Raine 1992). 
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2 Black hole thermodynamics 



Within purely classical framework, black holes in any diffeomorphism covariant theory of gravity 
(where the field equations directly follow from the diffeomorphism covariant Lagrangian) and in 
general relativity mathematically resembles some aspects of classical thermodynamic systems (Wald 
1984, 1994, 2001, Kiefer 1998, Brown 1995 and references therein). In early seventies, a series of 
influential works (Bekenstein 1972, 1972a, 1973, 1975, Israel 1976, Bardeen, Carter & Hawking 1977, 
see also Bekenstein 1980 for a review) revealed the idea that classical black holes in general relativity, 
obey certain laws which bear remarkable analogy to the ordinary laws of classical thermodynamics. 
Such analogy between black hole mechanics and ordinary thermodynamics ('the generalized Second 
Law', as it is customarily called) leads to the idea of the 'surface gravity' of black hole, 1 k, which 
can be obtained by computing the norm of the gradient of the norms of the Killing fields evaluated 
at the stationary black hole horizon, and is constant on the horizon (analogous to the constancy of 
temperature T on a body in thermal equilibrium - the Zeroth Law of classical thermodynamics). 
Also, k=0 can not be reached by performing finite number of operations (analogous to the 'weak 
version' of the third law of classical thermodynamics where temperature of a system cannot be made 
to reach at absolute zero, see discussions at Kiefer 1998). It was found by analogy via black hole 
uniqueness theorem (see. e.g. Heusler 1996 and references therein) that the role of entropy in classical 
thermodynamic system is played by a constant multiple of the surface area of a classical black hole. 



3 Hawking Radiation 

The resemblance between the laws of ordinary thermodynamics to those of black hole mechanics 
were, however, initially regarded as purely formal. This is because, the physical temperature of a 
black hole is absolute zero (see, e.g. Wald 2001). Hence physical relationship between the surface 
gravity of the black hole, and the temperature of a classical thermodynamic system can not be con- 
ceived. This further indicates that a classical black hole can never radiate. However, introduction of 
quantum effects might bring a radical change to the situation. In an epoch making paper published 
in 1975, Hawking (1975) used quantum field theoretic calculation on curved spacetime to show that 
the physical temperature and entropy of black hole does have finite non-zero value (see Page 2004 for 
an excellent review of black hole thermodynamics and Hawking radiation). A classical space time 
describing gravitational collapse leading to the Schwarzschild black hole was assumed to be the dy- 
namical back ground, and a linear quantum field, initially in it's vacuum state prior to the collapse, 
was considered to propagate against this background. The vacuum expectation value of the energy 
momentum tensor of this field turned out to be negative near the horizon. This phenomenon leads to 
the flux of negative energy into the hole. Such negative energy flux would decrease the mass of the 
black hole and would lead to the fact that the quantum state of the outgoing mode of the field would 
contain particles. 2 . The expected number of such particles would correspond to radiation from a 
perfect black body of finite size. Hence the spectrum of such radiation is thermal in nature, and the 
temperature of such radiation, the Hawking temperature T# from a Schwarzschild black hole, can be 
computed as 

Th ~ 87rk b GM BH ^ 

where c and kj, are the velocity of light in vacuum and the Boltzmann's constant, respectively. The 
semi classical description for Hawking radiation treats the gravitational field classically and the quan- 

1 The surface gravity may be defined as the acceleration measured by red-shifts of light rays passing close to the horizon 
(Heifer 2003) 

2 For a lucid description of the physical interpretation of Hawking radiation, see, e.g., Wald 1994, Kiefer 1998, Heifer 
2003. 
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tized radiation field satisfies the d'Alembert equation. At any time, black hole evaporation is an 
adiabatic process if the residual mass of the hole at that time is larger than the Planck mass. 

4 Towards an analogy of Hawking effect: The motivation 

Substituting the values of the fundamental constant in eq. (1), one can rewrite T H for a Schwarzschild 
black hole as : 

T "~ 6 - 2xl0 "(^)° A ' (2) 

It is evident from the above equation that for one solar mass black hole, the value of the Hawking 
temperature would be too small to be experimentally detected. A rough estimate shows that T# for 
stellar mass black holes would be around 10 7 times colder than the cosmic microwave background 
radiation. The situation for super massive black hole will be much more worse, as Thoc ^/ . Hence 
T H would be a measurable quantity for primordial black holes with very small size and mass, if such 
black holes really exist, and if instruments can be fabricated to detect them. The lower bound of mass 
for such black holes may be estimated analytically. The time-scale T over which the mass of the black 
hole changes significantly due to the Hawking's process may be obtained as (Heifer 2003) 

r~(^)V (3 ) 

As the above time scale is a measure of the lifetime of the hole itself, the lower bound for a primordial 
hole may be obtained by setting T equal to the present age of the Universe. Hence the lower bound 
for the mass of the primordial black holes comes out to be around 10 15 gm. The size of such a black 
hole would be of the order of 10~ 13 cm and the corresponding Th would be about 10 llo K, which is 
comparable, as we will see in §10, with the macroscopic fluid temperature of the freely falling matter 
(spherically symmetric accretion) onto an one solar mass isolated Schwarzschild black hole. How- 
ever, present day instrumental technique is far from efficient to detect these primordial black holes 
with such an extremely small dimension, if such holes exist at all in first place. Hence, the observa- 
tional manifestation of Hawking radiation seems to be practically impossible. 

On the other hand, due to the infinite redshift caused by the event horizon, the initial configuration 
of the emergent Hawking Quanta is supposed to possess trans-Planckian frequencies and it's wave 
lengths are beyond the Planck scale; thus low energy effective theories cannot self consistently deal 
with the Hawking radiation (Parentani 2002). Also, the nature of the fundamental degrees of freedom 
and the physics of such ultra short distance is yet to be well understood. Hence the fundamental is- 
sues like the statistical meaning of the black hole entropy, or the exact physical origin of the out going 
mode of the quantum field, remains unresolved (Wald 2001). 

Perhaps the above mentioned issues served as the principal motivations to launch a theory, analo- 
gous to the Hawking's one, whose effects would be possible to comprehend through relatively more 
perceivable physical systems. The theory of analogue Hawking radiation opens up the possibility to 
experimentally verify some basic features of black hole physics by creating the sonic horizons in the 
laboratory. A number of works have been carried out to formulate the condensed matter or optical 
analogs of event horizons 3 It is also expected that analogue Hawking radiation may find important 

3 Literatures on study of analogue systems in condensed matter or optics are quite large in numbers. Condensed matter 
or optical analogue systems deserve the right to be discussed as separate review articles on its own. In this article, we, 
by no means, are able to provide the complete list of references for theoretical or experimental works on such systems. 
However, to have an idea on the analogue effects in condensed matter or optical systems, readers are refereed to the book 
by Novello, Visser & Volovik 2002 for review, and to some of the representative papers like Leonhard 2002, 2003, Garay, 
Anglin, Cirac, & Zoller 2000, 2001, Jacobson & Volovik 1998, Volovik 1999, 2000, 2001, Brevik & Haines 2002, Schiitzhold, 
Giinter & Gerhard 2002, Schiitzhold & Unruh 2002, de Lorenci, Klippert & Obukhov 2003, Reznik 2000, Novello, Perez 
Bergliaffa, Salim, DeLorenci & Klippert 2003. 
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uses in the fields of FRW cosmology (Barcelo, Liberati & Visser 2003), inflationary models, quantum 
gravity and sub-Planckian models of string theory (Parentani 2002). 

For space limitation, in this article, we will, however, mainly describe the formalism behind the clas- 
sical analogue systems. By 'classical analogue systems' we refer to the examples where the analogue 
effects are studied in classical fluids, and not in quantum fluids. In the following section, we narrate 
the basic features of a classical analogue system. Hereafter, we shall use Tah to denote the analogue 
Hawking temperature, and T# to denote the the actual Hawking temperature as defined in eq. (1). 
We shall also use the words 'analogue', 'acoustic' and 'sonic' synonymously in describing the hori- 
zons or black holes. Also the phrases 'analogue Hawking radiation/effect/temperature' should be 
taken as identical in meaning with the phrase 'analogue radiation/ effect/ temperature'. A system 
manifesting the effects of analogue radiation, will be termed as analogue system. 



5 Theory of analogue radiation 

In a pioneering work, Unruh (1981) showed that a classical system, relatively more clearly perceivable 
than a quantum black hole system, does exist, which resembles the black hole as far as the quantum 
thermal radiation is concerned. The behaviour of a linear quantum field in a classical gravitational 
field was simulated by the propagation of acoustic disturbance in a convergent fluid flow. In such a 
system, it is possible to study the effect of the reaction of the quantum field on it's own mode of prop- 
agation and to contemplate the experimental investigation of the thermal emission mechanism. If one 
considers the equation of motion for a transonic barotropic irrotational fluid, it can be shown (Unruh 
1981) that the scaler field representing the acoustic perturbation (i.e, the propagation of sound wave) 
satisfies a differential equation which is analogous to the equation of a massless scaler field propa- 
gating in a metric. Such a metric closely resembles the Schwarzschild metric near the horizon. Thus 
acoustic propagation through a supersonic fluid forms an analogue of event horizon, as the 'acoustic 
horizon' at the transonic point. The behaviour of the normal modes near the acoustic horizon indi- 
cates that the acoustic wave with a quasi-thermal spectrum will be emitted from the acoustic horizon 
and the temperature of such acoustic emission may be calculated as: 



h 

Tah = 



1 du\ 
a, dn 



(4) 

acoustic horizon 



Where a s is the sound speed, tti is the component of the dynamical flow velocity normal to the 
acoustic horizon, and ^ represents the normal derivative. Equation (4) has clear resemblance with 
eq. (1) and hence Tah is designated as analogue Hawking Temperature and such quasi-thermal ra- 
diation form acoustic (analogue) black hole is known as the Analogue Hawking radiation. Note that 
the sound speed in Unruh's original treatment (the above equation) was assumed to be constant in 
space, i.e., equation of state used was isothermal. 

Unruh's work was followed by other important papers (Jacobson 1991, Unruh, 1995, Visser, 1998, 
Jacobson 1999, Bilic 1999). A generalized treatment for classical analogue radiation for Newtonian 
fluid may be obtained in Visser (1998). Visser (1998) showed that for barotropic, inviscid, hydrody- 
namic fluid, the equation of motion for the velocity potential ip describing an acoustic disturbance 
is identical to the d'Alembertian equation of motion for a minimally coupled massless scaler field 
propagating in a (3+1) Lorenzian geometry as: 

= ^La M (v=£ <r = o (5) 
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The acoustic propagation is described by the following metric, which is algebraically dependent on 
the density and dynamical as well as sonic velocity of the flow : 




(6) 



The above acoustic metric for a point sink, is conformally related to the Painleve-Gullstrand-Lemaitre 
form of Scwarzschild metric. The conformal factor may be neglected in connection to the calculation 
of the analogue Hawking temperature. The surface gravity k for above mentioned acoustic metric 
can be calculated as: 



From eq (7), it is straight forward to calculate the expression for analogue temperature as 



The above equation shows that one needs to know the exact location of the acoustic horizon, the ex- 
act values of the dynamical and sound velocities, and their space gradients on the acoustic horizon; 
to determine the analogue Hawking temperature of a classical analogue system. Note that the only 
difference between eq (8) and eq. (4) is that, in Unruh's original expression, positional constancy of 
sound speed was assumed, whereas such assumption has not been made in eq. (8). 
For analogue systems discussed above, the fluid particles are coupled to the flat metric of Mankowski's 
space (because the governing equation for fluid dynamics in the above treatment is completely New- 
tonian), whereas the sound wave propagating through the non-relativistic fluid is coupled to the 
curved pseudo-Riemannian metric. Phonons (quanta of acoustic perturbations) are the null geodesies, 
which generate the null surface, i.e., the acoustic horizon. Introduction of viscosity may destroy the 
Lorenzian invariance and hence the acoustic analogue is best observed in a vorticity free completely 
dissipationless fluid. That is why, the Fermi superfluids and the Bose-Einstein condensates are ideal 
to simulate the analogue effects. The most important issue emerging out of the above discussions is 
the following: 

Even if the governing equation for fluid flow is completely non-relativistic (Newtonian), the acoustic 
fluctuations embedded into it are described by a curved pseudo-Riemannian geometry. This infor- 
mation is useful to portray the immense importance of the study of the acoustic black holes, i.e. the 
black hole analogues, or simply, the analogue systems. 

In summary, analogue (acoustic) black holes (or systems) are fluid-dynamic analogue of general rel- 
ativistic black holes. Analogue black holes possess analogue (acoustic) event horizons at local tran- 
sonic points. Analogue black holes (alternatively, analogue systems) emit analogue Hawking radi- 
ation, the temperature of which is termed as analogue Hawking temperature, and it's expression is 
formulated using Newtonian description of fluid flow. Black hole analogues are important to study 
because it may be possible to create them experimentally in laboratories to study some properties 
of the black hole event horizon, and to study the experimental manifestation of Hawking radiation. 
In the subsequent sections, we will describe how accretion processes onto astrophysical black holes 
may constitute a concrete example of analogue system naturally found in the Universe. We will 
also demonstrate that accreting black holes are unique as analogue systems because they are the 
only system found in nature so far, where the analogue temperature may exceed the actual Hawking 
temperature. 




(8) 
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6 Transonic black hole accretion in astrophysics 



Gravitational capture of surrounding fluid by massive astrophysical objects is known as accretion. 
There remains a major difference between black hole accretion and accretion onto other cosmic ob- 
jects including neutron stars and white dwarfs. For celestial bodies other than black holes, infall of 
matter terminates either by a direct collision with the hard surface of the accretor or with the outer 
boundary of the magneto-sphere, resulting the luminosity through energy release from the surface. 
Whereas for black hole accretion, matter ultimately dives through the event horizon from where ra- 
diation is prohibited to escape according to the rule of classical general relativity and the emergence 
of luminosity occurs on the way towards the black hole event horizon. The efficiency of accretion 
process may be thought as a measure of the fractional conversion of gravitational binding energy of 
matter to the emergent radiation and is considerably high for black hole accretion compared to ac- 
cretion onto any other astrophysical objects. Hence accretion onto classical astrophysical black holes 
has been recognized as a fundamental phenomena of increasing importance in relativistic and high 
energy astrophysics. The extraction of gravitational energy from the black hole accretion is believed 
to power the energy generation mechanism of X-ray binaries and of the most luminous objects of the 
Universe, the Quasars and active galactic nuclei (Frank, King & Raine 1992). The black hole accretion 
is, thus, the most appealing way through which the all pervading power of gravity is explicitly man- 
ifested. 

If the instantaneous dynamical velocity and local acoustic velocity of the accreting fluid, moving 
along a space curve parameterized by r, are u(r) and a s (r) respectively, then the local Mach number 
M(r) of the fluid can be defined as M(r) = ^pj- The flow will be locally subsonic or supersonic ac- 
cording to M(r) < 1 or > 1, i.e., according to u(r) < a s (r) or u(r) > a s (r). The flow is transonic if at 
any moment it crosses M = 1. This happens when a subsonic to supersonic or supersonic to subsonic 
transition takes place either continuously or discontinuously. The point(s) where such crossing takes 
place continuously is (are) called sonic point(s), and where such crossing takes place discontinuously 
are called shocks or discontinuities. In order to satisfy the inner boundary conditions imposed by the 
event horizon, accretion onto black holes exhibit transonic properties in general. 
Investigation of accretion processes onto celestial objects was initiated by Hoyle & Littleton (1939) 
by computing the rate at which pressure-less matter would be captured by a moving star. Subse- 
quently, theory of stationary, spherically symmetric and transonic hydrodynamic accretion of adi- 
abatic fluid onto a gravitating astrophysical object at rest was formulated in a seminal paper by 
Bondi (1952) using purely Newtonian potential and by including the pressure effect of the accret- 
ing material. Later on, Michel (1972) discussed fully general relativistic polytropic accretion on to a 
Schwarzschild black hole by formulating the governing equations for steady spherical flow of per- 
fect fluid in Schwarzschild metric. Following Michel's relativistic generalization of Bondi's treat- 
ment, Begelman (1978) and Moncrief (1980) discussed some aspects of the critical points of the flow 
for such an accretion. Spherical accretion and wind in general relativity have also been considered 
using equations of state other than the polytropic one and by incorporating various radiative pro- 
cesses (Shapiro, 1973a,b, Blumenthal & Mathews 1976, Brinkmann 1980). Malec (1999) provided the 
solution for general relativistic spherical accretion with and without back reaction, and showed that 
relativistic effects enhance mass accretion when back reaction is neglected. The exact values of dy- 
namical and thermodynamic accretion variables on the sonic surface, and at extreme close vicinity of 
the black hole event horizons, have recently been calculated using complete general relativistic (Das 
2002) as well as pseudo general relativistic (Das & Sarkar 2001) treatments. 

For flow of matter with non-zero angular momentum density, spherical symmetry is broken and ac- 
cretion phenomena is studied employing axisymmetric configuration. Accreting matter is thrown 
into circular orbits around the central accretor, leading to the formation of the accretion discs around 
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the galactic and extra-galactic black holes. For certain values of the intrinsic angular momentum den- 
sity of accreting material, the number of sonic point, unlike spherical accretion, may exceed one, and 
accretion is called 'multi-transonic'. Study of such multi-transonicity was initiated by Abramowicz 
& Zurek (1981). Subsequently, multi-transonic accretion disc in general relativity has been studied in 
a number of works (Fukue 1987, Chakrabarti 1990, Kafatos & Yang 1994, Yang & Kafatos 1995, Pariev 
1996, Peitz & Appl 1997, Lu, Yu, Yuan & Young 1997, Das 2004, Barai, Das & Wiita 2004, hereafter 
BDW). All the above works, except BDW, usually deal with low angular momentum sub-Keplerian 
prograde flow. BDW studied the retrograde flows as well and showed that a higher angular mo- 
mentum (as high as Keplerian) retrograde flow can also produce multi-transonicity. Sub-Keplerian 
weakly rotating flows are exhibited in various physical situations, such as detached binary systems 
fed by accretion from OB stellar winds (Illarionov & Sunyaev 1975; Liang & Nolan 1984), semi- 
detached low-mass non-magnetic binaries (Bisikalo et al. 1998), and super-massive BHs fed by ac- 
cretion from slowly rotating central stellar clusters (Illarionov 1988; Ho 1999 and references therein). 
Even for a standard Keplerian accretion disk, turbulence may produce such low angular momentum 
flow (e.g., Igumenshchev & Abramowicz 1999, and references therein). 



7 Motivation to study the analogue behaviour of transonic black hole ac- 
cretion 

Since the publication of the seminal paper by Bondi in 1952, the transonic behaviour of accreting 
fluid onto compact astrophysical objects has been extensively studied in the astrophysics community, 
and the pioneering work by Unruh in 1981 initiated a substantial number of works in the theory 
of analogue Hawking radiation with diverse fields of application stated in §4. It is surprising that 
no attempt was made to bridge these two categories of research, astrophysical black hole accretion 
and the theory of analogue Hawking radiation, by providing a self-consistent study of analogue 
Hawking radiation for real astrophysical fluid flows, i.e., by establishing the fact that accreting black 
holes can be considered as a natural example of analogue system. Since both the theory of transonic 
astrophysical accretion and the theory of analogue Hawking radiation stem from almost exactly the 
same physics, the propagation of a transonic fluid with acoustic disturbances embedded into it, it 
is important to study analogue Hawking radiation for transonic accretion onto astrophysical black 
holes and to compute Tah for such accretion. 

In the following sections, we will describe the details of the transonic accretion and will show how 
the accreting black hole system can be considered as a classical analogue system. Hereafter, we 
define gravitational radius r g = 2G ^ BH , where Mbh is the mass of the black hole, G is the universal 
gravitational constant and c is the velocity of light. The radial distances and velocities are scaled in 
units of r g and c respectively and all other derived quantities are scaled accordingly; G = c = Mbh = 
1 is used. 



8 Newtonian and semi-Newtonian spherical accretion 

In this section we will follow the methodology developed in Das & Sarkar 2001 to describe the sim- 
plest possible example of transonic black hole accretion. If the gravitational potential is defined as 
then the equation of motion for spherically accreting matter onto the black hole can be written as: 

du du 1 dp d& 
dt dr p dr dr 

where symbols have their usual meaning. The first term of eqn. (9) is the Eulerian time derivative 
of the dynamical velocity at a given r, the second term is the 'advective' term, the third term is 
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the momentum deposition due to pressure gradient and the final term is due to the gravitational 
acceleration. Note that $ may be purely Newtonian, or may represent any of the following semi- 
Newtonian pseudo-Schwarzschild black hole potentials: 



2(r-l)' 



-1 + 



$9 = 



1 

2r 



3 / 1 

1 + 121 — 

2r \2r 



1 

1 - - 

r 



* A = hn(l-- 
2 V r 



(10) 



$i is proposed by Paczyhski and Wiita (1980), $2 by Nowak and Wagoner (1991), $3 and $4 by 
Artemova, Bjornsson & Novikov (1996). Introduction of such potentials allows one to investigate 
the complicated physical processes taking place in accretion in a semi-Newtonian framework by 
avoiding pure general relativistic calculations. Most of the features of space-time around a com- 
pact object are retained and some crucial properties of the analogous relativistic accretion solutions 
could be reproduced with high accuracy. Hence, those potentials might be designated as 'pseudo- 
Schwarzschild' potentials, see Das & Sarkar 2001 and Das 2002 for details. 
The continuity equation can be written as 



dp Id/ 2 \ 



(ii) 



For a polytropic equation of state, the steady state solutions of eqn. (9) and eqn. (11) provides: 
1) Conservation of specific energy £ of the flow: 



u 



+ 



2 7-I 



(12) 



where / y(=c p /c v ) is the adiabatic index of the accreting material. 
2) Conservation of Baryon number (or accretion rate M): 

M = 4npur 2 
1 1 

Defining the entropy accretion rate E = Mji- 1 K 1- 1 , we have: 

1 

S = Ana] _1 ur 2 
The space gradient of dynamical velocity is obtained as: 



du 
dr 



r 

u 



(13) 



(14) 



(15 



u 



where $ is equal to d^/dr. Since the flow is assumed to be smooth everywhere, if the denominator 
of eqn. (15) vanishes at any radial distance r, the numerator must also vanish there to maintain the 
continuity of the flow. One thus calculates the sonic point quantities as 



u 



(r=r h ) 



M(r=r h ) 







2 


(r=r h ) 



(16) 



where r h is the sonic point, a family of which forms the spherical acoustic horizon. The location of 
the acoustic horizon, i.e., r h , can be obtained by algebraically solving the following equation: 



(r=r h ) 







(17) 



(r=r h ) 



8 



Ur J ( r=rh ) 

tion: 



at its corresponding sonic point can be obtained by solving the following quadratic equa- 



(1 + 7) ( 



f du\- 
dr J 



+ 



(r=r h ) 



2.829 (7- 1) ^ 









(r=r h ) 





du\ 
dr J 



(r=r h ) 



+ (27-1) 



(r=r h ) 



= (18) 



(r=r h ) 



From the above discussion it is evident that for completely Newtonian flow = — as well as for 
semi-Newtonian flow, one obtains the exact location of the transonic point, which is the location of 
the acoustic horizon because of the stationarity approximation. One can also calculate the exact value 
of [u, a s ,du/dr, da s /dr] at the acoustic horizon. From eq. (8), it is obvious that all the required quan- 
tities to calculate the value of analogue Hawking temperature may be obtained from first principle 
from the calculations presented above, if the gravitational potential is purely Newtonian in nature. 
Hence, Newtonian spherical accretion onto astrophysical black holes is the example of analogue sys- 
tem where the exact value of Tah can he calculated using only two measurable parameters, namely, £ 
and 7. Also to be noted that the above formulation was carried out using position dependent sound 
speed. Unruh's original calculation of Tah (see eq. (4)) assumes the positional constancy of sound 
speed, which may not always be a real physical situation. This was modified to include the position 
dependent sound speed by Visser (1998) for flat space and by Bilic (1999) for curved space. However, 
none of these works provides the exact calculation of and acoustic horizon related quantities from 
first principles, and hence the numerical value of Tah in any existing literature (even for flat space, 
let alone curved space) was obtained using approximate, order of magnitude, calculations only. 
The most important point to mention in this context is that, the accreting astrophysical black holes are 
the only real physical candidates for which both the horizons, actual (electromagnetic) and analogue 
(sonic), may exist together simultaneously. Hence our application of analogue Hawking radiation to 
the theory of transonic astrophysical accretion may be useful to compare some properties of these 
two kind of horizons, by theoretically studying and comparing the behaviour of the same flow close 
to these horizons. Spherically accreting black hole system is, therefore, not only a concrete example 
of analogue system naturally found in the Universe, it is also most important among all the existing 
analogue models, because study of black hole accretion as black hole analogue may bridge the gap 
between two apparently different school of thoughts, the theory of transonic astrophysical accretion, 
and the theory of analogue Hawking radiation, sharing intrinsically the same physical origin. 
If $ represents any of the pseudo-Schwarzschild potentials defined by eq. (10), the analogue Hawk- 
ing temperature can still be determined from the knowledge of eq. (16 - 18). One only needs to 
modify eq. (8) by incorporating such potentials instead of the purely Newtonian one, while calculat- 
ing the acoustic metric and other related quantities. It is not difficult to show that (Das & Dasgupta 
2004, in preparation) for purely Newtonian as well as for semi-Newtonian black hole potentials, Tah 
can never exceed Th, which is not true for general relativistic accretion, as we will see in §10. 



9 General relativistic black hole accretion 



In this section, we will describe fully general relativistic astrophysical accretion, both spherically sym- 
metric and axisymmetric, onto astrophysical black holes. The mass of the accreting fluid is assumed 
to be much less compared to Mbh (which is usually reality for astrophysical black hole accretion), so 
that the gravity field is controlled essentially by Mbh only. Calculations presented in the following 
sections are mostly adopted from Das 2004, Das 2004a and BDW. 
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9.1 Mono-transonic spherical accretion 



Accretion is 9 and <f> symmetric and possesses only radial inflow velocity. Stationary solutions are 
considered. Accretion is governed by the radial part of the general relativistic time independent 
Euler and continuity equations in Schwarzschild metric. The conserved specific flow energy £ (the 
relativistic analogue of Bernoulli's constant) along each stream line reads £ = hu t (Anderson 1989) 
where h and are the specific enthalpy and the four velocity which can be re-cast in terms of the 
radial three velocity u and the polytropic sound speed a s to obtain: 



7-1 



7 -(l + a2 



1 



u 



2 ' 



(19) 



We concentrate on positive Bernoulli constant solutions. The mass accretion rate M may be obtained 
by integrating the continuity equation: 



M = 4-n-pur 



1 



r (1 — u 2 ) ' 



(20) 



where p is the proper mass density. A polytropic equation of state, p = Kp~< , is assumed (this is 
common in astrophysics to describe relativistic accretion) to define E as: 



E = K f- 1 M = Anpur 



r (1 — v?) 



a 2 s (7 ~ 1) 
7-(l + af) 



(21) 



Simultaneous solution of eq. (19 - 21) provides the dynamical three velocity gradient at any radial 
distance r: 



du 
dr 



u (1 - u 2 ) [a 2 s (4r - 3) - 1] 



2r (r — 1) (u 2 



The sonic point conditions are obtained as: 



u 



(r=r h ) 



s|(r=r h ) 



, 477, - 3 

Substitution of u|( r=rh ) and a s |( r=rh ) into eq. (19) for r = r h provides: 



where 



Ti 



r\ + r^Tx + r h T 2 + T 3 = 0, 
~2£ 2 (2-37) + 9(7- 1) 



(22) 



(23) 



(24) 



4( 7 -l)(£2-i) 
£ 2 (3 7 - 2) 2 - 27( 7 -l) 2 



, r 3 



27 



64 {£ 2 - 1)' 



(25) 



32 (£ 2 - 1)( 7 - l) 2 

Solution of eq. (24) provides the location of the acoustic horizon in terms of only two accretion 
parameters {£ ,7}, which is the two parameter input set to study the flow. Eq. (24) can be solved 
completely analytically by employing the Cardano-Tartaglia-del Ferro technique. One defines: 

3T 2 -T? „ 9TiT 2 - 27T 3 - 2T? o o 
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n x = {/S2 + JSI + S?, 
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so that the three roots for come out to be: 



(26) 



*r h = -^--±(p--iV3to-). (27) 

However, note that not all V/j{i = 1, 2, 3} would be real for all {£, 7}. It is easy to show that if * > 0, 
only one root is real; if 4* = 0, all roots are real and at least two of them are identical; and if * < 0, 
all roots are real and distinct. Selection of the real physical (r h has to be greater than unity) roots 
requires the following discussion. 

Although one can analytically calculate and other variables at rh, there is no way that one can ana- 
lytically calculate the flow variables at any arbitrary r. One needs to integrate eq. (24) numerically to 
obtain the transonic profile of the flow for all range of r, starting from infinity and ending on to the 
actual event horizon. To do so, one must set the appropriate limits on {£, 7} to model the realistic sit- 
uations encountered in astrophysics. As £ is scaled in terms of the rest mass energy and includes the 
rest mass energy, £ < 1 corresponds to the negative energy accretion state where radiative extraction 
of rest mass energy from the fluid is required. For such extraction to be made possible, the accret- 
ing fluid has to possess viscosity or other dissipative mechanisms, which may violate the Lorenzian 
invariance. On the other hand, although almost any £ > 1 is mathematically allowed, large values 
of £ represents flows starting from infinity with extremely high thermal energy, and £ > 2 accre- 
tion represents enormously hot flow configurations at very large distance from the black hole, which 
are not properly conceivable in realistic astrophysical situations. Hence one sets 1 ^ £ ^ 2. Now, 
7 = 1 corresponds to isothermal accretion where accreting fluid remains optically thin. This is the 
physical lower limit for 7; 7 < 1 is not realistic in accretion astrophysics. On the other hand, 7 > 2 
is possible only for superdense matter with substantially large magnetic field (which requires the 
accreting material to be governed by general relativistic magneto-hydrodynamic equations, dealing 
with which is beyond the scope of this article) and direction dependent anisotropic pressure. One 
thus sets 1 ;$ 7 ^ 2 as well, so {£, 7}has the boundaries 1 ^ {£, 7} ;$ 2. However, one should note 
that the most preferred values of 7 for realistic black hole accretion ranges from 4/3 to 5/3 (Frank et. 
al. 1992). 

Coming back to the solution for rh, one finds that for the preferred range of {£, 7}, one always obtains 
^ < 0. Hence the roots are always real and three real unequal roots can be computed as: 

\ h = 2v^cos (§) - T -±?r h = 2^ cos (e + 2?[) _ ^ 

\ h = 2v^cos (®±^) - £ (28) 

One finds that for all 1 < {£, 7} < 2, \^ becomes negative. It is observed that {V/j, 3 ^} > 1 for most 
values of the astrophysically tuned {£ , 7}. However, it is also found that \h does not allow steady 
physical flows to pass through it; either u, or a s , or both, becomes superluminal before the flow 
reaches the actual event horizon, or the Mach number profile shows intrinsic fluctuations for r < r h . 
This information is obtained by numerically integrating the complete flow profile passing through 
\h- Hence it turns out that one needs to concentrate only on for realistic astrophysical black hole 
accretion. Both large £ and large 7 enhance the thermal energy of the flow so that the accreting fluid 
acquires the large radial velocity to overcome a s only in the close vicinity of the black hole . Hence 
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anti-correlates with {£, 7}. To obtain (dv/dr) and (da s /dr) on the acoustic horizon, L' Hospital's rule 
is applied to eq. (22) to have: 

= $ 12 - $123, 

r=r h 



dr 



~-if - 1 I * 12 * 123 ) (29) 
dr L =ri \ JW h - 3 2 2" 



r=r h 



where 
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6r h (r h -l) 

1 Vi^3 ' 2 477,- 3 

$3 = ^"H fr-g (7 " I) 2 " r h (10 7 2 - 19 7 + 9) 
(4r h -3)2 V > 

+ (6 7 2 - 11 7 + 3)], * 4 = '^V^-l^"^ (30) 
9.2 Multi-transonic disc accretion 

In this section, we follow the treatment provided by Das 2004 and BDW because only in those works 
a consistent methodology has been developed which can handle the transonic solutions upto the ex- 
treme close vicinity of the black hole event horizon (down to a distance upto around 1.001r fl ). The 
Boyer-Lindquist coordinates with signature — h ++ is used, and an azimuthally Lorentz boosted 
orthonormal tetrad basis co-rotating with the accreting fluid. A is defined to be the specific angu- 
lar momentum of the flow and any gravo-magneto-viscous non-alignment between A, and the Kerr 
parameter a, is neglected. Stationary axisymmetric solution of the following equations is considered: 

VS^^O, {pv") =0, (31) 

where Of^", v^, and p are the energy momentum tensor, the four velocity and the rest mass density 
of the accreting fluid, respectively. The semicolon denotes the covariant derivative. The 'stationar- 
ity' condition implies the vanishing of the temporal derivative of any scalar field in the disc or the 
vanishing of the Lie Derivative of any vector or tensor field along the Killing vector Jj. This implies 
that |j2 = if II is a scalar field, and Cq/qiII. = if II represents a vector or a tensor field. The 

'axisymmetry' property is endowed with the space like Killing field ^ J^J M . 

The radial momentum balance condition may be obtained from (v^v v + g^ u ) 9^" = 0. Exact so- 
lutions of these conservation equations require knowledge of the accretion geometry and the in- 
troduction of a suitable equation of state. Polytropic accretion is considered for which p = Kp" 1 , 
where K and 7 are the monotonic and continuous functions of the specific entropy density, and the 
constant adiabatic index of the flow, respectively. The specific proper flow enthalpy is taken to be 
h = (7 — 1) {7 — (l + a 2 )} 1 , where a s is the polytropic sound speed defined as a s = {dp/de) l J 2 = 
*i (T(r), 7) = *2 (p, P, 7); here T(r) is the local flow temperature, e is the mass-energy density, and 
^2} are known functions. The subscript S indicates that the derivative is taken at constant spe- 
cific entropy. Hence, a 2 (r) = jnT(r) /{pmu) = @ 2 T(r), where = ^/(/im/;)] 1 / 2 , p is the mean 
molecular weight, run is the mass of the hydrogen atom, and k is Boltzmann's constant. 
The temporal component of the first part of eq. (31) leads to the conservation of specific flow energy 
£ . The Kerr metric in the equatorial plane of the black hole may be written as (e.g., Novikov & Thorne 
1973) 

ds 2 = g^dx^dx" = - T -^-dt 2 + — (# - oodtf + T —dr 2 + dz 2 , (32) 
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where A = r — 2r + a , A = r 4 + r 2 a + 2ra , and u> = 2ar/A. The angular velocity SI is 



n = - 



4a A(4a'-r'A) 
+ *9tt) _ — a 



4Aa 



The normalization relation v^v^ = — 1, along with the above value of $7, provides: 

_ i 

Ar 2 A 



(33) 



vt 



(1 - u 2 ) {A 2 - 8XarA + X 2 r 2 (4a 2 - r 2 A)} 



(34) 



where u is the radial three velocity in the co-rotating fluid frame. Hence the conserved specific energy 
(which includes the rest mass energy) can be re-written as: 



(7-1) 



[7 _ ( 1 + 02 T)] 



Ar 2 A 



(1 - n 2 ) {A 2 - 8XarA + A 2 r 2 (4a 2 - r 2 A)} 



(35) 



For a suitably chosen disc geometry (see BDW for details about the disk structure in Kerr metric, 
also see Das 2004, for the details of relativistic accretion in a different disc geometry in Schwarzschild 
metric) The mass and entropy accretion rate comes out to be: 



M = AirpQT^ Mr 2 
H = A L \p2■Kr 2 u^ 



2A (7- 1)6 2 T 



7 (i - g 2 tm 2 ) {7 - (i + e 2 r)} ip 

l(l±l\ 

2 \-f-l) 



(1-u 2 )^ 

The velocity gradient can be calculated as: 



gj (7 ~ 1) 
7 { 7 -(l + a 2 )} 



(36) 
(37) 



du 
dr 



2a\ 
(7+1) 


r-1 1 2 vtcrx 
A ' r Aip 


X 
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u 2a| 


1 - 


u 2 v t rJ 


(l-M 2 ) (7+1)(1-« 2 )m 


2ip 



where 



a = 2X 2 v t -a 2 , X = + 



IdA X dn (%^ + A%0 



A dr (1 - QX) dr (g^ + A^) 
The sonic point conditions can be obtained as: 



(38) 



(39) 



a s \( 



r=r h ) 



u 2 (7 + 1) j) 
2tp — U 2 VtCF 



_, 1 

2 



J (r=r h ) 



l l(: 



r=r h ) 



2r (r - 1) + 4A 



(r=r h ) 



(40) 



For any value of [£, X, 7, a], substitution of the values of u|( r=rh ) and a s |( r=rh ) in terms of r^ in the 
expression for £ (see eq. (35)), provides a polynomial in r^, the solution of which determines the loca- 
tion of the sonic point(s) r^. After lengthy algebraic manipulations, the following quadratic equation 
is formed, which can be solved to obtain (du / dr) ( r =r h ) : 



a 



du\ 



+ P 



/ du\ 



dr J (r=Ph) \drj {r=rh) 



+ C = o, 



(41) 



where the coefficients are 
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Q: 



{1 + U 2 
,2^2 



{i-u 2 y 7 + 1 

2 



, P = — — + r 6, C = -T5; 



7 + 1 



#4 = (52 



a 2 (1 - 6 2 ) A 
5l " «(1- U 2)' 52 

, #5 = 



2 m; t <5 3 



U 2 Vt<J 

" 2^ ' 
3n 2 -l 



1 2A 2 a 
63 = - + -, 

v t <J tp 



u ' 1 — u 2 
5 e = 



_ h ^ (7 ~ 1 ~ aj) 
u{l-u 2 ) 1-S 2 a 2 (l-u 2 ) 

(7 - 1 - a 2 ) x . fahxvt 



2a 2 



+ 



r - 1 2 crt; t x 



ar 2X 1 2 (r - 1)' 

T3 = ^' T4= A" 



A 2 



2(1-52)' 
(4A 2 vt - a 2 ) tp - v t a 2 
aip 

2_ _ v t adx 
r 2 4ip dr ' 



T5 



7 + 1 

T6 = ■ 



' « 2r 4 - { (7 - 1 - a 2 ) ti + ^a 2 r 3 } 

— (7 - 1 - a s J + a s T 3 



Idx 
2 dr' 



( 7 + l)(l-u 2 



(42) 



Note that all the above quantities are evaluated at r^. 

Putting the Kerr parameter equals to zero in the corresponding equations above, one can obtain the 
value of sonic point and [u,a s ,du/dr,da s /dr] at the sonic point, for accretion onto Schwarzschild 
black holes. 

Adopting a different disc geometry, the corresponding quantities in Schwarzschild metric can also 

directly be computed as the following (Das 2004): 

1) The sonic point is obtained by solving the following equation: 



where 



and 



4 + A 2 (1 - r h ) 



V (rh, A) 



l-*(r h ,A) 
7 + 1 



rh ~ 1 _ I" r h (7 - 1) 
7-f?(rh,A) 

I + (r h , A) 



= 



(43) 



* (r h , A) = 

2) The sonic velocities can be found as: 

/ 9 

«l(r=p h ) : 
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2r h -l 
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2r h (r h -l) 



l-2r h 



2r(r h -l) + 2[rl+^(l-r h )\ 
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' 1 
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2r 3 -A 2 


2r-l 1 


_r :1 +A 2 (l-7)_ 


2r(r-l) 


l-2r , A 2 -3r 2 
[ 2r(r-l) 1 2[r s +X I (l-r)\ J 



(44) 



(r=r h ) 



3) The velocity gradient can be obtained by solving the following equation 

"7 -(1 + 4) 



2(2 7 -3a 2 ) fdu\ 2 

(7 + DK-1) 2 UA + 4e(rh ' A) 



ut 



du\ 
dr)h 
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7 + 1 



a^(r h ,A) 



2£(r h ,A) 



7-(l + «I) 
7 + 1 



2r fe -l 

r fe (r h - 1) r 3 + A 2 (1 - r ft ) 



20r% - Ylr\ - 2A 2 (3r fe - 2) 
'5r£-4r 3 - A 2 (3^-4^ + 1) 



(45) 



where 



e(r h ,A) 



1 - 2r 



+ 



A 2 - 3r 2 



2r (r - 1) 2 [r 3 + A 2 (1 - r)] 



In the above two equations, and ah indicates the value of u and a s on the acoustic horizon (on 
r = rh), and the subscript h in general indicates that the quantities are measured on the acoustic 
horizon. The relation between the dynamical velocity gradient and acoustic velocity gradient at any 
point (including at the sonic point) can be obtained as: 



da s (r) 
dr 



7 + 



T {7-(l + a2)} 



du 



where: 



f(r,A) = 



1 - 2r 



+ 



u (u 2 — 1) dr 



A 2 - 3r 2 



+ f(r,A) 



(46) 



(47) 



2r (r - 1) ' 2 (r 3 - A 2 r + A 2 ) 

Note that unlike spherical accretion, one finds three sonic points in disc accretion for some values of 
[S , A, 7] for accretion onto the Schwarzschild as well as the Kerr black hole, among which the largest 
and the smallest values correspond to the X type outer, r D , and inner, r, u sonic points, respectively. 
The O type middle sonic point, r m , which is unphysical in the sense that no steady transonic solution 
passes through it, lies between r^ and r D . If the accretion through r D can be perturbed in a way so that 
it produces an amount of entropy exactly equal to [£(rj) — £(r )], supersonic flow through r Q can join 
the subsonic flow through rj by developing a standing shock. The exact location of such a shock, as 
well as the details of the post-shock flow, may be obtained by formulating and solving the general 
relativistic Rankine-Hugoniot condition for the flow geometry described above. Note that as a flow 
through ^ does not connect the black hole event horizon with infinity (such a flow folds back onto 
itself, see, e.g., Fig. 1 of BDW), solutions through rj do not have independent physical existence, and 
can only be accessed if the supersonic flow through r Q undergoes a shock, so that it generates extra 
entropy, becomes subsonic, and produces the physical segment of the solution through rj. 



10 Analogue Hawking temperature for general relativistic accretion 
10.1 Relativistic sonic geometry 

As mentioned in §8, it is not quite difficult to calculate Tah for Newtonian, and also, perhaps, for 
semi-Newtonian black hole accretion. For general relativistic flow, however, non-trivial effort has to 
be made to accomplish such goals. This is because, treatments in Unruh (1981) and Visser (1998), 
and in related works, are based on the distinct fact that the background fluid is governed by New- 
tonian spacetime, whereas for general relativistic accretion presented in §9.1 & 9.2 in this paper, the 
background fluid is characterized by purely Schwarzschild metric. Hence a relativistic version of 
acoustic geometry is necessary to calculate Tah for this purpose. The work by Bilic (1999) is the first 
(and perhaps, the only one in the literature as far as our knowledge is concerned) which calculates 
the generalized expression for acoustic analogue of surface gravity for propagation of relativistic 
transonic fluid. The relativistic acoustic metric comes out to be (Moncrief 1980, Bilic 1999): 



^ v ha* 



(l - afj u^] (48) 
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where h is the relativistic enthalpy density and N is the particle number density, 
surface gravity may be calculated as: 



The corresponding 



1 — a 2 Q dn 



(u 



(49) 



acoustic horizon 

Where is the stationary Killing field and d/dnis the normal derivative. 
10.2 Analogue temperature for spherical accretion 

For spherically symmetric general relativistic flow onto Schwarzschild black holes described in §9.1, 
one can evaluate the exact value of the Killing fields and Killing vectors to calculate the surface 
gravity for that geometry. The analogue Hawking temperature for such geometry comes out to be 
(Das 2004a): 



Tah 



he 3 



A-Kk b GM BH 



rj (r h - 0.75) 



d 
dr 



u) 



(50) 



r=r h 



where the values of rh, (du/dr)h and (da s /dr)h are obtained using the system of units and scaling 
used in this article. 

It is evident from the above formula that the exact value of Tah can be analytically calculated from 
the results obtained in §9.1. While eq. (27) provides the location of the acoustic horizon (rh), the 



value of 



dr 



(a — u) is obtained from eq. (29 - 30) as a function of £ and 7, both of which are 



r=r h 



real, physical, measurable quantities. Note again, that, since and other quantities appearing in 
eq. (50) are analytically calculated as a function of {£, 7}, eq. (50) provides an exact analytical value 
of the general relativistic analogue Hawking temperature for all possible solutions of an spherically 
accreting astrophysical black hole system, something which has never been done in the literature 
before. If \/Ar h - 3(1/2 - 1/3>4)(<J>12 - $123) > 1/ one always obtains (da s /dr < du/dr) h from eq. (29), 
which indicates the presence of the acoustic white holes at r^. This inequality holds good for certain 
astrophysically relevant range of {£, 7}; see following discussions. 

For a particular value of {£, 7}, one can define the quantity r to be the ratio of Tah and Th as: 



r = TahITh- 



(51) 



This ratio r comes out to be independent of the mass of the black hole, which enables us to compare 
the properties of two kind of horizons (actual and acoustic) for an spherically accreting black hole 
with any mass, starting from the primordial holes to the super massive black holes at galactic cen- 
ters. 

Note, however, that the analogue has been applied to describe the classical perturbation of the fluid 
in terms of a field satisfying the wave equation on a curved effective geometry. Main motivation of 
the methodology described in this section is not a rigorous demonstration of how the phonon field 
generated in this system could be made quantized. To accomplish that task, one can show that the 
effective action for the acoustic perturbation is equivalent to the field theoretical action in curved 
space, and the corresponding commutation and dispersion relations (see, e.g., Unruh & Schiitzhold 
2003) may directly follow from there. 

In the accompanying figure (taken from Das 2004a), the complete {£,7}space is classified accord- 
ing to the value of r. For accretion with {£, 7}taken from the lightly shaded regions marked by 
< Tah < Th, the Hawking temperature dominates over its analogue counterpart. For high 7, high 
£ flow ('hot' accretion), analogue Hawking radiation becomes the dominant process compared to the 
actual Hawking radiation (Tah > Th). For low £ and intermediate/high 7, acoustic white holes 
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Figure 1: Parameter space classification for {£, 7} for four different sets of values of Tah- 

appear (white region marked by WH, where da s /dr < du/dr at r h ). For {£, 7}belonging to this re- 
gion, one obtains outflow (outgoing solutions with Mach number increasing with increase of r) only. 
This has also been verified by obtaining the complete flow profile by integrating eq. (24). The dark 
shaded region (lower right corner of the figure) represents {£, 7} for which comes out to be physi- 
cal (rh > 1) but $f 23 becomes negative, hence (du/dr)h and (da s /dr)h are not real and Tah becomes 
imaginary. Note that both Tah > T H and white hole regions are obtained even for 4/3 < 7 < 5/3, 
which is the range of values of 7 for most realistic flows of matter around astrophysical black holes. 
Hence the domination of the analogue Hawking temperature over the actual Hawking temperature 
and the emergence of analogue white holes are real astrophysical phenomena. 

Note one important point that the accreting spherical black holes in general relativity is the only ana- 
logue system found till date, where the analogue temperature may be higher than the actual Hawking 
temperature. Newtonian or semi-Newtonian accreting system does not show this behaviour. Also 
note that both Th and Tah comes out to be quite less compared to the macroscopic classical fluid 
temperature of accreting matter 4 . Hence black hole accretion system may not be a good candidate to 
allow any observational test for detecting the analogue radiation. 

10.3 Analogue effects in accretion disc 

Calculation of Tah for disc accretion (as described in §9.2) is also possible. One needs to calculate the 
acoustic metric and the Killing vectors for axisymmetric space time. Then one needs to use those re- 
sults to calculate the surface gravity and analogue temperature. Initial attempts (Das & Bilic, 2004, in 
preparation) shows that although the geometry would be different (axisymmetric) than the spherical 
case, Tah will still be a function of [u, a s , du/dr, da s /dr] at the sonic point, with different functional 
form as that of eq. (50). Hence for disc accretion also, the exact value of Tah can directly be calculated 
from the results obtained in §9.2. However, as we discussed before, for axisymmetric accretion, more 
than one sonic point is formed, hence the question now may arise that which sonic points are to be 
taken, to form the acoustic horizon. The answer seems to be quite straight forward. Firstly, all the 
middle sonic points would be excluded because no stable solutions pass through that. Secondly, as 
accretion can pass through the inner sonic point only if a shock forms (see §9.2), one should exclude 
the inner sonic points as well. This is because, shock production is dissipative and may violate the 
Lorenzian invariance. Hence, if the accretion is multi-transonic, which is the situation for certain 
values of [£,\,j,a], only the flow passing through the outer sonic point will resemble an analogue 

4 The value of the flow temperature for non-radiative spherical accretion is found to be of the order of 10 10 — 10 11 °K . 
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system. For other values of [£, A, 7, a], disc accretion is mono-transonic and hence there will not be 
any confusion about the choice of the sonic points for producing the acoustic horizon. 

However, although an observer at positive infinity will see only one acoustic black hole, there may 
be one alternative possibility for an observer situated in between the two sonic points. An observer in 
between two sonic points (the inner sonic point and the outer sonic point) may perhaps see a white 
hole and a black hole or two white holes (see Barcelo et. al. 2004 for discussion about somewhat 
similar situation). Detail study of this issue is beyond the scope of this article and will be presented 
elsewhere. 
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